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1. Introduction 



Nuclear reactions govern major aspects of the chemical evolution of galaxies and stars 
(Fowler, 1984, Morel et al., 1999). The analytic study of the reaction rates and reaction 
probability integrals was undertaken by Critchfield (1972), Anderson et al. (1994), Haubold 
and Mathai (1998), and Chaudhry (1999). A proper understanding of the nuclear reactions 
that are going on in hot cosmic plasma, and those in the laboratories as well, requires a 
sound theory of nuclear-reaction dynamics (Clayton, 1983, Bergstroem et al., 1999). The 
rate of unlike reacting nuclei i and j in the case of nonrelativistic nuclear reactions taking 
place in nondegenerate environment is expressed as (Clayton, 1983, Lang, 1999) 

where rii and rij denote the particle number densities of the reacting nuclei i and j, /i — 

^ — — is the reduced mass of the reacting nuclei, T is the temperature, k is the Boltzmann 

rrii + rrij 

constant, (t{E) is the cross-section for the reaction under consideration, and v — {j^^ is 
the relative velocity (for reactions between like nuclei i — j, fiiTij has to be replaced by n^). 
Thus (Clayton, 1983, Lang, 1999) 

/•CO 

Tij = TiiHjX := niUjiav) := J a{E)v{E)i){E)dE (1.2) 

is the definition of the reaction probability integral A, that is, the probability per unit 
time that two nuclei, confined to a unit volume, will react with each other. The reaction 
probability is written in the significant form {av) to indicate that it is an appropriate average 
of the product of the reaction cross section and relative velocity of the interacting nuclei. If 
the reacting mixture is in thermal equilibrium this quantity depends only on the temperature. 

For nonresonant nuclear reactions between nuclei of charges Zi and Zj at low energies 

(below the Coulomb barrier), the reaction cross-section has the form (Clayton, 1983, Rowley 
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and Merchant, 1991, Lang, 1999) 



a{E) 



E 



(1.3) 



where 



m = (f ) 




(1.4) 



is the Sommerfeld parameter, h is Planck's quantum of action, and e is the quantum of 
electric charge. Eq. (1.3) allows the extrapolation of measured reaction cross sections down 
to astrophysical energies by introducing the S-factor. 

It is to be noted that S{E), a residual function of energy, represents intrinsically nuclear 
parts of the probabihty for the occurrence of a nuclear reaction (Clayton, 1983). It is often 
found to be constant or a slowly varying function over a limited energy range when the 
interaction energy of the pair of nuclei is not nearly equal to an energy at which the two 
nuclei resonate in a quasi-stationary state. 

The normalized energy distribution, representing the isotropic velocity distribution of the 
reacting nuclei, is given by (Clayton, 1983, Lang, 1999) 



(for reaction rates with anisotropic distributions see Imshennik, 1990). The substitution of 
ip{E)dE in (1.2) yields 




(1.5) 





where ^ = ( f 



1/2 



h 



,2 



. If the residual function S{E) is considered to be constant, 5*0, the 



corresponding reaction probability integral is given by 




1/2 



(1.7) 
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In effect, the reaction probability is governed by the average of the Gamow penetration factor 
over the Maxwell- Boltzmann distribution. Anderson et al. (1994) considered for S{E) the 
Maclaurin series expansion 

S{E) - S{0) + -^E + -^-^E (1.8) 

and solved (1.6) in terms of G- and i^-functions (Mathai and Saxena, 1973, 1978). 

It is to be noted that the residual function S{E) may not admit the representation (1.8) 
if the thermonuclear fusion plasma is not in the thermodynamic equilibrium, that is, the 
case in which there is a depletion or cut-off of the high energy tail of the Maxwell- Boltzman 
distribution in (1.5) (Lapenta and Quarati, 1993). Some of these cases have been considered 
in Anderson et al. (1994) and Haubold and Mathai (1998). In this paper we consider 
representations for the residual function given by 

S^{E) = So5{E-Eo), (1.9) 



S2{E) = So ) , (1.10) 



/ E 



a-l 



SsiE) = So[^^J H{E-E,), (1.11) 

Ey-^ 

—J exp{C{E-E^)), (1.13) 



S^E) - So[—J H(E^-E), (1.12) 



where 



1, if E> El, 
H(E-Ei):^i (1.14) 
^ 0, if E < Ei 



is the unit step function and 



5(E-Ei):=^{H(E-Ei)) (1.15) 



is the Dirac delta function (Jennings and Karataglidid, 1998). 
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We discuss the analytic representations of the corresponding reaction probabihty inte- 
grals. The integral corresponding to (1.9) may be called instantaneous reaction probability 
integral. 

2. The Astrophysical ThermonucleEir Functions 

The theoretical and experimental verification of nuclear cross-sections leads to the deriva- 
tion of the closed-form representation of the thermonuclear reaction rates. These rates are 
expressed in terms of the four astrophysical thermonuclear functions (given in the notation 
chosen by Anderson et al., 1994) 

h{z,v):= / y-exp(-y-— )dy, (2.1) 

h{z,d.y) /%"exp(-y-^)%, (2.2) 
Jo ^ 

h{z,t,v) ■= / y^exp(-y- ^==)dy, (2.3) 

h{z,6,b,u):= y^exp{-y-by'-—)dy, (2.4) 
Jo A/y 



„2 



where y = = ^ relates E or v, respectively, to the mean thermal velocity and z — 
(2^^) ^ — 27r (2^) ^ OiZiZj, were the velocity of light was introduced to make 

the dimension more apparent and to show the dependence on Sommerfeld's fine structure 
constant a. The closed-form representations of these integral functions, in terms of G- 
and H-functions, asymptotic values and, numerical results are discussed in Anderson et al. 
(1994) and Chaudhry (1999). Accounts about the developments of Meijer's G- function and 
of Fox's H-function are given in Mathai and Saxena (1973) and Mathai and Saxena (1978), 
respectively. Originally the investigations on these generalized hypergeometric functions were 
confined to theoretical results such as their analytical properties, integral representations, 
and asymptotic expansions, and then to the study of symmetric Fourier kernels, the solution 
of certain functional equations and other mathematical topics. Later, both the G-function 



and the H-function have been also used in the fields of statistical and astrophysical sciences, 
including extensive studies of the thermonuclear functions in eqs. (2.1)-(2.4) (Mathai, 1993). 



3. The Extended Gamma Functions 



The study of the astrophysical thermonuclear functions led to the development of a new 

class of special functions (Cliaudhry and Zubair, 1998). In particular, we define the extended 
gamma functions by 

r{a,x;b]f5) := H e-^e'*-^'^^ dt, (3.1) 

Jx 

and 

7(a,a;;6;/3) := r e-^e-'-^''^ dt. (3.2) 
Jq 

It is to be noted that the functions (3.1) and (3.2) are special cases of the general class of 
extended gamma functions introduced in Chaudhry and Zubair (1998). 



In fact we have 



pc+ioo 



ib,x) 



((0,1), (a, (3)) 
:= — / r(s)r(a + /3s, x)b-'ds, 



and 



'y{a,x;b;l3) = 70,2 



{b,x) 



((0,1), (a,/?)) 

]_ pc+ioo 

:= — / 7(s)r(a + ps, x)b~'ds. 
These functions satisfy the decomposition formula (Chaudhry and Zubair, 1998). 



^{a,x;b,P) + r{a,x;b;P)^H^g 



{b,x) 



(0,1), {a, (3)) 



(3.3) 



(3.4) 



(3.5) 



The astrophysical thermonuclear functions (2.1) - (2.4) are special cases of the extended 
gamma functions (3.1) and (3.2) (Chaudhry, 1999) 

h{z, ^, ^) = e* t ^ ) i-tr-'r + 1, t; z; 



" i-hY / 1\ 



'r=0 



(3.6) 
(3.7) 
(3.8) 

(3.9) 



It is straightforward to note that the transformation theorem 
1 



r(Q;, X] b; P) 



P' 



-.2,0 
0,2 



((o,^), («,1)) 



a; > 0,6 > 0,/3 > 0), (3.10) 



for the extended gamma function reveals as a special case that 



X a 





0-2,0 
-"0,2 


nj 


n+1,0 
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0,n+l 







(0,1), (a,n) 



(27r)(^-")/^v^ 



5 ) ■ ■ ■ ) 

1 2 n-l 

U, , , . . . , , Oi 

n n n 



(3.11) 



The relation (3.11) yields directly the closed form representation (Anderson et al., 1994) 

h{z,v) = v(v+l,Q-z^- 



'z" 




4 


0,^,1 + ^^ 



(3.12) 



The asymptotic representation of the extended gamma function V{ol^ 0; h\ /?) for small and 
large values of h are given as follows (Chaudhry and Zubair, 1998) 

T{pL) + -^r l^-^j + o(6), for small 6, a 0, 

for small 6, a = 0, 



r(a,0;6;/3) 



(3.13) 



In 6, 



and 



r(a,0;6;/3) 



X exp 



1 / 27r/3 



1/2 



^(2a+/3)/2(l+/3) ^(2a-l)/2(/3+l) 



^1 + ^)/3/(l+/3) 5l/(l+/3) 



for large 6. 



(3.14) 



The representations (3.13) and (3.14) are useful in finding the asymptotic values of the 
thermonuclear reaction probability integrals. 



4. Thermonuclecir Reaction Probability Integrals 



In this section we study the closed form representations of the thermonuclear reaction 
probability integrals corresponding to the various forms (1.9) - (1.13) of the residual function 
S{E). 

(4.1) Instantaneous Residual Function, S{E) — SoS{E — Eq). The substitution of this 
value of S{E) in (1.6) yields an elementary representation 

of the corresponding thermonuclear probability integral. 

(4.b) Power Type Residual Function, S{E) = So{E / Eo)""-^ . The substitution of the 
power type residual function in (1.6) yields 

/ o \ 1/2 



/;(^/^o)-exp(-A__L),^. (4.2) 



The transformation t = E / kT m. {A.2) yields 

which is available from Haubold and Mathai (1998) and can be written in terms of the 
extended gamma function, i.e.. 

,,.m-^f-r-'rLo,4,;iV (4.4) 



8 



In view of the decomposition relation (3.5) we can simplify (4.4) to obtain 



1/2 



A2= — 



kT 



that can further be simplified to give 



2,0 
0,2 



A, = — 



\ 'JO 



an I ^/Eq \ Eq 



(0,1), («,2) 



62 1 ■ 

ttt; 0, -, q; 
kT '2' 



(4.5) 



(4.6) 



(4.c) Power Type Delayed Residual Function, S{E) = Sq{E / Eq^-^ H {E - Ei). The 
substitution of the power type delayed residual function in (1.6) gives 



1/2 



A. = — 



5*0 



The transformation t = E/kT in (4.7) yields 



/;(W-exp(-A__L),^. (4.7) 



'/'^(kry-^^ / b_ 

^[EoJ JE./kT ''''^l ,/kT 



A, 



//TT 



kTViJ 



] dt, 



or 



(4.8) 



(4.9) 



The asymptotic representations of A3 for small and large values of 6/ can be found from 
(3.13) and (3.14). In fact for small value of h/\/¥T and a^Q 

2r(-a)6n 



A. 



'Eq \ Eq 



kT 



(4.10) 



(4.d) Power Type Cut-off Residual Function, S{E) = Sq{E/ Eq^-^H^E^ - E). The 
substitution of the power type cut-off residual function in (1.6) gives 



1/2 



A4 



j^E/E^y-^e-Kv^- 



E 



dE 



\lin J {kTy/^ Jo ^ ' kT ./E, 

that can similarly be simplified in terms of the extended gamma function to give 



(4.11) 



A4 



'kT' 



) ^/Eq \Eq 



2 / El 6 1 

7 «, 



kT' 2 



(4.12) 
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(4.e) Exponential Type Residual Function, S{E) = So{E / E^Y'^ e^^{-C {E - ^i)), 
{C > — The substitution of the exponential type residual function in (1.6) yields 



The transformation 



in (4.13) yields 



t= (c + -^)e, (c>--^ 

* kTj ' V kT 



(4.14) 

The integral in (4.14) is solvable in terms of extended gamma function to give 
The equation (4.15) can be simplified further in terms of the H-function to give 



5. Conclusion 



kT 



CkT+1 



(0,1), (a, 2) 



(4.16) 



Nuclear reactions govern major aspects of the chemical evolution of galaxies and stars. 
The analytic study of the reaction rates and reaction probability integrals is important in 
astrophysics. The development of the new class of extended gamma functions by Chaudhry 
and Zubair has facilitated in convenient notations for these analytic representations. We 
have considered various forms of the residual function occurring in the reaction probability 
integrals and have written them analytically in terms of the extended gamma functions. 
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